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1. Introduction 

IceCube is a cubic-kilometer neutrino detector at the Geographic South Pole capable to detect 
neutrinos in the energy regime above ^ lOGeV. It consists of 86 strings in the Antarctic glacial 
ice, each equipped with 60 photosensors deployed at depths between 1450m and 2450m below 
the surface [1]. The detection of high-energy neutrino events is based on the measurement of 
Cherenkov photons, emitted by charged particles produced by neutrino interactions in the ice. 
IceCube recorded the data used for this work from 2008 to 2011 in the detector configurations 
IC40, IC59 and /C79, where the numbers denote the amounts of deployed strings. 

Astrophysical neutrinos are well suitable messengers for studying production and accelera¬ 
tion mechanisms of cosmic rays. Recently, IceCube discovered a diffuse all-flavor astrophysical 
neutrino flux [2]. However, the properties of its astrophysical sources have not been identified 
yet. Besides single point source and catalogue searches [3], an angular correlation analysis was 
performed on IceCube data from 2008 to 2011 to search for sources too weak to be detected in¬ 
dividually [4]. In this work, the results of this analysis are reinterpreted with a new method to 
account for realistic source populations and the observed astrophysical neutrino flux. 

2. Method 

2.1 Limits used for astrophysical interpretation 

In this section, the key properties of the angular correlation analysis are summarized. It is 
based on spherical harmonics expansions of muon neutrino arrival direction skymaps. The signal 
expectation has been benchmarked according to different signal hypotheses, characterized by three 
quantities: The number of sources A^sou; the mean number of measured neutrinos per source at 
the horizon p {source strength)', and the astrophysical spectral index 7 . Assuming a constant flux 
from each source, a source at the horizon is assigned a number of neutrinos following a Poisson 
distribution with mean p; with changing declination, the mean number of neutrinos per source is 
reduced according to the lower detector acceptance. The signal part of each hypothesis includes 
Asou isotropically distributed sources in the Northern Sky with such Poisson distributed numbers 
of neutrinos. To match the total number of 108310 neutrinos in agreement with the number of 
events in the used IceCube data sample, a background part for the respective hypothesis is added. 
It consists of atmospheric neutrinos, sampled from experimental data. A test statistic (TS) value 
is calculated from the spherical harmonics expansion coefficients of each skymap. It denotes how 
significantly the specific skymap is distinguishable from pure atmospheric background. By repeat¬ 
ing this process for many skymaps of each specific signal hypothesis, a distribution of TS values is 
obtained for each signal hypothesis. The shift of this distribution with respect to the distribution of 
the TS of the purely atmospheric background in units of the standard deviation of the distribution 
of the TS of the background is called signalness Z. Finally, a Feldmann-Cousin approach [5] is 
used to compute experimental limits (90 % C.L.) in terms of Asou, M and 7 from the TS value of 
the experimental skymap and the distributions of the TS of the simulated hypotheses. 

2.2 Motivation and conversion 

While the angular correlation analysis described in Section 2.1 tests signal hypotheses for fixed 


2 



Neutrino angular correlations in IceCube 


Michael S chimp 



Figure 1: Signalness per source against 
source strength for astrophysical energy spec¬ 
tra and legend: best-fit power law 

exponent 



log(M) 

Figure 2: Sketch of the source count distribution 
with powers j8i and P 2 adopted from 
Fermi-LAT 


values of Asou and /i, corresponding to a fixed flux per source at Earth, realistic source populations 
contain multiple sources of variable strength. Thus, to obtain limits on realistic source distribu¬ 
tions, two possibilities exist. First, every source population of interest can be simulated and its 
distribution of the TS investigated; second, the already obtained limits can be converted to limits 
on realistic source populations. To describe a realistic source population, a so called source count 
distribution is used which describes the number of sources per (infinitesimal) source-strength 
interval. The source count distribution is a function of the signal parameter /i, used in the angular 
correlation analysis, which simplifies the conversion of the limits from the analysis to limits on 
realistic source distributions as shown in the following. It should be noted that while this work 
focuses on the limits, also other signalness-related results from the angular correlation analysis, 
like the sensitivity or the discovery potential, can be converted in the same way. To perform the 
conversion of the limits, two quantities have to be equated: the signalness of the limit Ziim from 
the angular correlation analysis (Ziim = 1.38 [4]), and the signalness of the source population of 
interest, given on the right hand side of Equation (2.1): 



0 


^^^(jU) is the source count distribution and is the signalness per source which is known 

from [6] to follow a power law with an exponent compatible with about 2 as shown in Figure 1. 
Thus, the signalness of the source population is the integral over the source strength ji of all sources 
per ji interval times their respective signalness per source. Solving Equation (2.1) for parameters 
of an assumed source count distribution yields the converted limits on these parameters. 

2.3 Exemplary source count distribution parametrization adopted from Fermi-LAT 

This work focuses on extragalactic neutrino sources. The gamma-ray telescope Fermi-LAT 
has detected extragalactic high-energy photon sources at high galactic latitudes [7] with a fitted 
energy spectrum of E~'^-^ which are also candidates for high-energy neutrino sources [8]. Thus, for 
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the neutrino source count distribution in this work, the parametrization of the gamma ray source 
count distribution from Fermi-LAT is adopted. It should be noted that this parametrization is, 
despite its physical motivation, only an exemplary parametrization to demonstrate the performance 
of this method. 

The Fermi-LAT source count distribution from [7] is the photon flux dependent number density 
^^^(5) where Asou is the number of sources and S is the photon flux per source at Earth. It is 
parametrized by a broken power law: 

cWsou _ \ A5-^S5>Sb 

dS 5<5b 


Pi = 2.49 zb 0.12 and P 2 = 1.58 zb 0.08 are the powers of the source count distribution after and 
before the break, respectively. A = (16.46 zb 0.80)10“^'^cm^deg“^s is a simple scale factor for 
and 5b = (6.60 zb 0.91)10“^cm“^s“^ is the particle flux at the break of the source count 
distribution. For the purpose of this work, we express this parametrization in terms of the (neutrino) 
source strength as defined in Section 2.1 (the source strength at the break is denoted /ib) instead of 
the photon flux S (and a breaking flux 5b) and in terms of a dimensionless scale factor B instead of 
A. All parameters we use are dimensionless in contrast to the replaced parameters. Furthermore, a 
cutoff is introduced by setting the source count distribution to zero for all source strengths above 
a maximum /imax because otherwise contributions from sources with diverging strength would be 
included in the signalness calculation (Equation (2.1)). It is fixed to /imax = 20/ib corresponding 
to the brightest source in the Fermi-LAT sample [7] which has a flux of 5max ~ 205b. The best-fit 
values for the powers of the Fermi-LAT source count distribution, j 8 i = 2.49 and P 2 = 1.58, are 
adopted for neutrinos (except for the generalization in Section 3.4). This results in: 


dAsou 

dju 


0, /i > 20/ib 

20/ib > jtt > Atb 


Figure 2 shows an illustration of the used parametrization. 


(2.3) 


3. Application and results 

The interpretation of the IceCube results proceeds according to the following steps: In Section 
3.1, Equation (2.1) is solved for the parameters B and /ib of the source count distribution described 
in Section 2.3 using the limit signalness from the angular correlation analysis. In Section 3.2, the 
parameters of the source count distribution are constrained to reproduce the expected numbers of 
neutrinos from the observed astrophysical neutrino flux. In Section 3.3, the Fermi-LAT observation 
is expressed in the parameters B and /ib as a function of a universal neutrino-to-photon ratio 
from all sources measured by Fermi-LAT. Values for e^jy corresponding to the measurements of 
Sections 3.1 and 3.2 are determined. As a last step, in Section 3.4, the dependencies of the results 
for Eyjy on Pi and P 2 are investigated. For the applications in Sections 3.1, 3.2 and 3.3, the powers 
of the source count distribution, j 8 i and P 2 are fixed to the Fermi-LAT best-fit values from Section 
2.3. The neutrino energy spectral indices 7 = 2.0 and 7 = 2.5 are tested, motivated by IceCube 
measurements of the astrophysical neutrino flux [ 2 ]. 
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Figure 3: Dashed lines: IceCube limits converted to 5(jUb); colored solid lines: (5, /ib) reproducing 
the observed astrophysical neutrino flux (s. Section 3.2); triangles: (5,jUb) at limit for universal 
Ey/y, diamonds: prediction for Eyjy reproducing the observed astrophysical neutrino flux purely 
originating from high-latitude sources measured by Fermi-LAT; asterisks: (5,jUb) for e^jy = 1; 
black line: allowed (5,/ib) for universal £y/y 

3.1 Limit conversion 

The solution of Equation (2.1) with from Equation (2.3) is a function 5(jUb). It represents 
the upper limit from the angular correlation analysis [4] on B for each value of /ib and is shown as 
a dashed line for each spectral index in Figure 3. The negative slopes of the exclusions lines are ex¬ 
plained by the increased signalness for larger /Tmax ^ Mb (s. Equation (2.1)) due to a corresponding 
larger non-zero integration range. This increased signalness is compensated by lower values for B, 
causing the negative slope. 

3.2 Constraints from the observed astrophysical neutrino flux 

The astrophysical neutrino flux, observed by IceCube in 2014 [2], can also be represented 
in terms of B and Mb- The number of measured signal neutrinos expected from the source count 
distribution fz(scd) and the number expected from the observed flux fi(astroph.) are equated: 

n(scd) = J dll 1 „(astroph.) = — (3.1) 

where IC denotes the IceCube detector configuration, the livetime of IceCube for the operation 
of each detector configuration, is the effective area of each configuration and ^ the observed 
differential astrophysical neutrino flux by IceCube. / is a factor ^ 1 to obtain the declination- 
averaged number of measured neutrinos for a value of ji. The pairs (5, Mb) solving this equation 
for the different considered neutrino energy spectra are the values for B and Mb that correspond to 
the number of neutrinos expected from the observed astrophysical neutrino flux and are shown as 
colored lines in Figure 3. They are different since Az(astroph.) is energy spectrum dependent due 
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7 ey/y prediction 8 ^/^ limit ratio 
092 OTS 0.29 

2.5 41.4 30.6 1.35 

Table 1: Results for universal Ey/y, 2nd column: predictions assuming the observed neutrino 

flux; 3rd column: Eyjy (upper) limits; 4th column: ratio between e^jy prediction and e^jy limit 


to the correlation of the flux normalization and the energy spectrum. Their intersections with the 
dashed limit lines (of the respective energy spectrum) separate the allowed regions (below) from 
the excluded regions (above) of parameter values constituting the observed flux. 


3.3 The Fermi-LAT best-fit value and the neutrino-to-photon ratio 


Assuming a universal neutrino-to-photon ratio at Earth e^jy from all sources measured in the 
Fermi-LAT high-latitude survey, one can express the gamma-ray source count distribution mea¬ 
sured by Fermi-LAT in terms of B, /ib (i.e. properties of a neutrino flux) and e^jy. First, the values 
for B and /ib corresponding to an equal flux of neutrinos and photons at Earth {e^/y = 1) in the en¬ 
ergy range of the used IceCube measurement (lOOGeV — lOOTeV), have to be determined. These 
particular values for B and /ib are called ^Permi and /ib,Fermi- Fermi-LAT measured the extragalactic 
photon flux in the energy range from 100MeV to lOOGeV. Thus, this particle flux is extrapolated 
into the energy range of IceCube using the assumed energy spectra with indices 7 = 2.0 and 7 = 2.5, 
yielding the values for ^Permi and /ib^Permi- The extrapolation from the energy range of Fermi-LAT 
into the energy range of IceCube yields different flux normalizations for the energy spectra and 
therefore different values for ^Permi and /ib^Permi- For oach value of 8 y/y, the source strength ji of 
each source in the population with B = ^Permi and jU}^ = Mb,Fermi has to be multiplied by £y/y. For 
the source count distribution parameters B and Mb. this leads to: 


B = £ 


i8i-l 

v/r 


•^Fe 


Mb — jy' Mb,Fermi 


(3.2) 


Thus, for each value of a certain pair (5, Mb) is obtained, resulting in the black line shown in 
Figure 3. The intersections of the (dashed) limit lines from Section 3.1 and the (solid) lines rep¬ 
resenting the observed astrophysical neutrino flux with the (black) line representing the neutrino- 
to-photon ratios 8 y/y directly yield upper limits on 8 y/y and predictions for 8 y/y for each neutrino 
energy spectrum in this simplified model. The values for the 8 y/y limit and the 8 y/y prediction are 
given in Table 1. 


3.4 Changing parameter space - varying source count distribution powers 

The procedure from Section 3.1 through 3.3 is repeated with varied powers of the source 
count distribution j 8 i and J 82 (defined as in Equation (2.3)). Figures 4a and 4b show the results 
in terms of the 8 y/y prediction as introduced in Section 3.2 and Figures 4c and 4d show the 8 y/y 
limit as introduced in Section 3.3 for the neutrino energy spectra with 7 = 2.0 and 7 = 2.5. Figures 
4e and 4f show the ratios between the 8 y/y prediction and the 8 y/y limit. Ratios larger than 1 
indicate that the astrophysical neutrino flux is excluded with 90 % C.L. to origin purely from source 
populations with the corresponding source count distributions with universal 8 y/y in the respective 
energy ranges denoted in Section 3.3. 
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Figure 4: (a), (b): predicted universal neutrino-to-photon ratios 8 ^/^ for source populations that 
correspond to the observed astrophysical neutrino flux for different powers jSi, j 82 (s. Equation 
(2.3)); (c), (d): converted limits on universal for different powers j 8 i, J 82 (s. Equation (2.3)); 
(e), (f): ratios between the e^jy prediction and the Eyjy limit; asterisk with error bars: j 8 i, J 82 and 
uncertainties from Fermi-LAT [7]; black line: color-coded quantity = 1 


4. Conclusions 

The universal neutrino-to-photon ratios Eyjy resulting in a number of neutrinos corresponding 
to the observed astrophysical neutrino flux by IceCube, called the Eyjy prediction, and the limit on 
set by the angular correlation analysis in [4] are summarized in Table 1. For y— 2.0, the ob¬ 
served astrophysical neutrino flux corresponds to a lower Eyjy prediction than the e^jy limit. Thus, 
for this energy spectrum, the sources from the Fermi-LAT high-latitude survey are not excluded 
to be the origin of the astrophysical neutrino flux under the stated assumptions. For y— 2.5, the 
opposite is the case: The Eyjy prediction exceeds the Eyjy limit and is thus excluded. 


7 
























































Neutrino angular correlations in IceCube 


Michael S chimp 


The results of the variations of j 8 i and J 82 (s. Equations (2.2) and (2.3)) reveal several insights: 
First, j 8 i plays a major role for both, the prediction and the limit since the source count 
distribution depends strongly on /3i for all source strengths ji. Second, the £y^y limit is almost 
independent of j 82 as shown in Figures 4c and 4d. As j 82 only affects the source count distribution 
for source strengths below the source strength of the break /ib (s. Equation (2.3)), this leads to 
the conclusion that the sources brighter than /ib are the signalness dominating sources. In contrast 
to that, for the Ey^y prediction corresponding to the astrophysical neutrino flux, there is a remarkable 
dependence on j 82 . The quantity used for the £y/y prediction is the number of neutrinos from the 
tested source count distribution fi(scd). As the integrand in the formula for fi(scd) (Equation (3.1)) 
is one power of ji lower compared to the integrand in the formula for the signalness (Equation 
( 2 . 1 )), fainter sources (i.e. small ji) affect the £y^y prediction while barely affecting the £y^y limit. 
A ratio between the £y^y prediction and the £y^y limit larger than 1 is excluded with 90 % C.L. 
Thus, for both examined energy ranges, the areas to the lower left from the black lines in Figures 
4e and 4f are excluded. For 7 = 2.5, where the hypothesis with the benchmark values for /3i and 
P 2 is excluded, their uncertainty intervals reach into the allowed region. 

Acknowledgments 

This work is supported by the Federal Ministry of Education and Research (BMBF), the 
Helmholtz Alliance for Astroparticle Physics (HAP) and the German Research Foundation 
(DFG). We thank Markus Ahlers, Jan Auffenberg and Leif Radel for the valuable discussions 
during the preparation of this manuscript. 

References 

11] A. Achterberg et al. (IceCube Collaboration), First Year Performance of The IceCube Neutrino 
Telescope, Astropart. Phys. 26 155 (2010) larXiv: astro-ph/0 60 4 450] 

12] M. G. Aartsen et al. (IceCube Collaboration), Evidence for Astrophysical Muon Neutrinos from the 
Northern Sky with IceCube, PRL (accepted) larXi v :1507.04005 [ as t ro-ph . HE ] ] 

13] M. G. Aartsen et al. (IceCube Collaboration), Searches for Extended and Point-like Neutrino Sources 
with Four Years of IceCube Data, ApJ 196 109 (2014) larXiv: 1406.6757 [ astro-ph. HE ] ] 

14] M. G. Aartsen et al. (IceCube Collaboration), Searches for small-scale anisotropies from neutrino 
point sources with three years of IceCube data, Astroparticle Physics 66 39 (2015) 

larXiv:1408.0634 [astro-ph.HE] ] 

15] G. J. Feldmann, R. D. Cousins, Unified approach to the classical statistical analysis of small signals, 
Phys. Rev. D 57 3873 (1998) larXiv:physics/9711021] 

16] M. Leuermann, Search for diffuse neutrino point sources using a multipole analysis in IceCube, 
Master’s thesis, RWTHAachen University (2014) 

17] A. A. Abdo et al., The Eermi-LAT high-latitude Survey: Source Count Distributions and the Origin of 
the Extragalactic Diffuse Background, ApJ 720 435 (2010) larXiv :1003.0895 

[astro-ph.CO]] 

[8] M. Petropoulou, S. Dimitrakoudis, P. Padovani, A. Mastichiadis, E. Resconi, Photohadronic origin of 
y-ray BL Lac emission: implications for IceCube neutrinos, Mon. Not. Roy. Astron. Soc. 448 2412 
(2015) [arXiv:1501.07115 [astro-ph . HE] ] 


8 




